Competing magnetic anisotropies in atomic-scale junctions 
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Using first-principles calculations, we study the magnetism of 5d transition-metal atomic junctions 
including structural relaxations and spin-orbit coupling. Upon stretching monatomic chains of W, 
Ir, and Pt suspended between two leads, we find the development of strong magnetism and large 
values of the magnetocrystalline anisotropy energy (MAE) of up to 30 meV per chain atom. We 
predict that switches of the easy magnetization axis of the nanocontacts upon elongation should 
be observable by ballistic anisotropic magnetoresistance measurements. Due to the different local 
symmetry, the contributions to the MAE of the central chain atoms and chain atoms in the vicinity 
of the leads can have opposite signs which reduces the total MAE. We demonstrate that this effect 
occurs independent of the chain length or geometry of the electrodes. 
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I. INTRODUCTION 

Fascinating insights into the formation of atomic 
chains consisting only of a few atoms suspended be- 
tween two electrodes have been obtained by transmission 
electron microscopy and mechanically controllable break 
junction techniques These experiments triggered the 
imagination of scientists to use such atomic-scale junc- 
tions as future electronic devices by exploiting their 
unique properties, e.g. ballistic electronic transport .^^"^ 
Due to their reduced dimensions, atomic chains have 
also been predicted to develop magnetic moments even 
for elements nonmagnetic in bulk such as Ir, Pt, or 
Pdi^i^ Recently, an indirect proof has been given that 
transition-metal chains in break junctions are in general 
magnetic.'^ Mastering the intriguing magnetic properties 
of suspended chains would enable spintronic applications 
based on the unique possibility to probe, control and 
switch the magnetic states by spin-polarized currents. ^iS 

The orientation of the magnetic moments is stabi- 
lized against thermal fluctuations by the magnetocrys- 
talline anisotropy energy (MAE) arising from spin-orbit 
coupling. Theoretical studies of this key quantity per- 
formed for idealized systems such as free-standing infinite 
monowires (MWs) and small clusters suggest giant val- 
ues for Ad- and Sd-transition-metals and switching of the 
easy axis upon stretching the wiresJ^^— Even the size 
of the magnetic moment itself can crucially depend on 
the magnetization direction, an effect coined as colossal 
magnetic anisotropy.— The ballistic conductance in such 
wires varies with the orientation of the magnetization di- 
rection with respect to the chain axis giving rise to ballis- 
tic anisotropic magnetoresistance (BAMR).^-^ While the 
predicted effects probably occur in real atomic-scale junc- 
tions, the theoretical studies have either focused on ide- 
alized systems or neglected varying interatomic distances 
in the chains.—"— However, for small suspended chains 
variations in the interatomic distance due to stretching 



the electrodes are unavoidable and can provide a unique 
way to control the properties of the system, not available 
in higher dimensions i^iiiiii 

Here, we go beyond such theoretical studies by ex- 
plicitly taking the contacts into account and by inves- 
tigating the development of magnetism in atomic-scale 
junctions upon increasing the distance between the elec- 
trodes. We perform first-principles calculations for W, Ir 
and Pt junctions including spin-orbit coupling and struc- 
tural relaxations and focus on atomic chains of three 
atoms suspended between two bulk- like bcc-(OOl) elec- 
trodes as shown in Fig. [1] These transition-metals are 
common tip materials in scanning tunneling microscopy, 
and Ir and Pt have been shown to form long atomic chains 
in break junction experiments.^ In particular we concen- 
trate on the influence of varying tip-to-tip separation on 
the magnetic properties of the junctions. We observe that 
relaxations of the chain atoms upon stretching are in ac- 
cordance to observed expermentally trends of chain for- 
mation for these elements.— We also find that large spin 
moments develop in the suspended trimers upon stretch- 
ing and on the central atom they reach similar values as 
in infinite monowires.—"— Interestingly, the W trimer be- 
haves as a single spin impurity as only the central atom 
develops a sizable magnetic moment. 

By including spin-orbit coupling (SOC), we obtain the 
MAE of the suspended chains and observe large values 
of 10—30 meV per atom. Surprisingly, we find that the 
contribution to the MAE from central chain atoms and 
atoms in contact with the electrodes can have opposite 
signs and are of similar giant magnitude of a few tens 
meV per atom. The magnitude of this effect, unknown 
in bulk or on surfaces may lead to a strong tendency to- 
wards non-coUinear magnetic order in atomic-scale junc- 
tions even when the exchange interactions are small or 
profoundly collinear,^-^ We prove the generality of this 
phenomenon by showing that it occurs for chains of differ- 
ent length with different geometry of the contacts. Fur- 
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FIG. 1: (color online) Upper panel: three-atoms chain sus- 
pended between two bcc(001)-contacts. L denotes the length 
of the unit cell consisting of 17 atoms, while da and dc give 
the distance between the first plane of the contact to the 
apex atom and from the apex atom to the central atom of the 
trimer, respectively. Lower panel: relaxation of the dc and da 
bonds in suspended trimers of W, Ir, and Pt upon increasing 
L. For comparison, the bond length for ideal stretching with 
fixed apex atom positions is given (dashed line). 



ther, the MAE may exhibit changes of sign upon stretch- 
ing, i.e. the easy magnetization direction switches from 
along the chain to perpendicular to it. We demonstrate 
that switching may be observable experimentally by mea- 
suring the ballistic conductance of the junctions. 



II. METHOD AND STRUCTURE 



netic ground state, while W trimer always converges to 
the antiferromagnetic solution. 

Spin-orbit coupling SOC was included into the calcu- 
lations in a non-perturbative way.^^ We define the mag- 
netic anisotropy energy (MAE) as the total energy differ- 
ence between configurations with a magnetization along 
the chain axis (z-direction), and perpendicular to it (r- 
direction) . The MAE was calculated employing the force 
theorem, and for several electrode stretchings, AL, we 
checked the force theorem results for the trimer 's MAE 
against self-consistent calculations, finding good qualita- 
tive agreement. Taking into account the magnitude of 
the obtained MAE we can safely neglect the magnetic 
dipolc-dipole contribution. 

In order to check the applicability of our ID geometry 
for the electrodes, which allows us to efficiently perform 
the demanding calculations with structural relaxations 
and including SOC, we have also considered a geometry 
with electrodes represented by surfaces, c.f. Fig. 2(a). 
Due to the immense computational cost, we have re- 
stricted our test to the most critical case of a Pt trimer 
and chose an elongation of AL — 2.4 a.u. For these 
calculations we used the film version of the FLEUR code. 
We performed a single super-cell calculation in which the 
trimer was suspended between two electrodes consisting 
of a four Pt atom square base deposited onto three lay- 
ers of fee Pt(OOl). This configuration will be referred to 
as "2D" in the following. Fig. 2. The ideal bulk values 
of Pt with an in-plane lattice constant of 5.23 a.u. were 
taken for the four contact atoms and the bulk-like leads. 
The positions of the trimer's atoms were set to their re- 
laxed values at AL = 2.4 a.u. obtained within the wire 
geometry of Fig. 1. The in-plane separation between the 
neighboring trimers was set to 15.7 a.u. We used a basis 
functions cut-off of fcmax = 3.7 a.u.~^ and 15 fc-points in 
the irreducible part of the 2D Brillouin zone, while 144 
fc-points in the whole Brillouin zone were used in order 
to determine the MAE with the force theorem. Along 
the trimer's axis only the F-point was used. 



We have performed first-principles calculations within 
the generalized gradient approximation (GGA) to the 
density functional theory (DFT). We employed the full- 
potential linearized augmented plane-wave method for 
one-dimensional (ID) systems (ID-FLAPW)r^ as imple- 
mented in the FLEUR codci^ Basis functions were ex- 
panded up to fcmax = 4.0 a.u.~^ and we have used 8 
A;-points in one half of the ID Brillouin zone. Consid- 
ering more fc-points does not change the values of the 
spin moments and magnetic anisotropy energies signifi- 
cantly. The contact structure has first been relaxed sep- 
arately in all three dimensions and then kept fixed for 
all contact separations L (Fig. 1). The apex atoms of 
the trimer have been relaxed along the chain direction 
{z) upon varying L, while the central atom is fixed by 
symmetry (Fig. [1}. Calculations have been performed 
for junctions stretched from the equilibrium up to AL — 
3 a.u. For Pt and Ir trimers we considered the ferromag- 



III. RELAXATIONS AND MAGNETISM 

Under the condition of reduced coordination of transi- 
tion metal atoms in a monoatomic chain suspended be- 
tween thicker leads, a situation realized e.g. in a break 
junction experiment, the distance d between the chain 
atoms becomes of utter importance. Upon pulling apart 
the leads during the chain elongation process^iii d is con- 
stantly changing which has dramatic consequences for the 
magnetic properties of suspended chains, in particular for 
heavy Ad and 5d transition-metalsi^i^iii By increasing 
the separation AL between the leads and accounting at 
the same time for the changes in the chain's interatomic 
distances — which was neglected in most of the previous 
ab initio studies of magnetism in suspended chains — we 
aim at mimicking a real break junction experiment. 

The relaxations of the trimers along the chain's sym- 
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FIG. 2: (color online) 2D configuration, and Bl and B2 cases 
of the Pt trimer considered for the separation of Pt electrodes 
of AL = 2.4 and 1.8 a.u., respectively. In 2D configuration 
the trimer is suspended between two electrodes consisting of 
a four Pt atom square base deposited onto three layers of fee 
Pt(OOl). Green color in Bl and B2 marks the Pt atoms of the 
contacts in which the spin moment was artificially quenched. 
In accordance to Fig. 1, red color marks the central atom and 
the blue color marks the apex atoms of suspended trimer. For 
more details see text. 



metry z-axis are shown in Fig. [l] where the bond lengths 
dc and da are presented as a function of contact separa- 
tion L. For W with an almost half-filled 5(i-shell and large 
bulk cohesive energy almost the entire increase in contact 
separation enters the dc bond length, while the da bond 
between the apex atom and the contact is nearly unaf- 
fected. In contrast, for Ir and Pt, some of the stretching 
transfers into an increased bond length at the contact. 
The apex atom of the Ir trimer has a stronger bond to 
the contact than to the central atom, which is reflected in 
the linear regime of dc{L) for AL > 1.5 a.u. For Pt, the 
dc bond is significantly stronger resulting in a remarkable 
stretching of the da bond even at large electrode separa- 
tions, which in a realistic break junction will eventually 
result in chain elongationJ^ Overall, our predicted trend 
is in accordance with experimental observations of in- 
creased probability for chain formation when going from 

W to Pt.2'22 

The magnetic properties of suspended 5c? TM atoms 
are very sensitive to the local environment. In most pre- 
vious studies it was common to model the contacts by 
semi-infinite surface electrodesJ^"— This restriction on 
the geometry of the contacts is quite strong in particular 
for 5d TM break junctions, as the surface geometry leads 
to suppressed spin moments of the atoms in suspended 
wire and contact atoms closest to it.^"* However, it is 
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FIG. 3: (color online) Local spin moments of the central (/ig) 
and apex (/ig) atoms (denoted as "c-atom" and "a-atom", 
respectively) of the suspended (a) W, (b) Ir and (c) Pt trimers 
as a function of the stretching of the contacts (without SOC). 
For comparison the values for an infinite free monowire (MW) 
are given by dashed lines. The interatomic distance in the 
infinite MW corresponds to the distance dc for a given AL. 
For Pt triangles down and triangles up stand for the values of 
the apex (filled) and central (open) spin moments in B2 and 
2D cases, respectively. For details see text. 



well-established that the shape of the contacts in a break 
junction is rather a thinning wire-like geometry, with the 
reduced coordination number of the contact atoms adja- 
cent to the chain which greatly enhances their tendency 
towards magnetism. In this respect our contact ge- 
ometry is probably closer to the real situation, although 
it might overestimate the tendency of the contacts to 
magnetism. In order to investigate this point quantita- 
tively, we have performed calculations in the setup shown 
in Fig. 2, which are discussed further below. 

The calculated spin moments p^s inside the atomic 
spheres of the suspended trimers upon increasing the dis- 
tance between the contacts AL are shown in Fig. [3] For 
W, the central atom reveals a sizable magnetic moment 
of 1.7/iB already at AL = which further increases upon 
stretching, while the moments of all other atoms are neg- 
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ligible. In this respect the W trimer presents a unique 
system of a single spin impurity between nonmagnetic 
leads. At small stretching, the spin moments of the Ir 
trimer are rather small with significantly larger moments 
of the central atom. Upon further stretching the apex 
and central moments rise rapidly, reaching large values of 
1.4/iB and respectively. Both, the spin moments 
of the central W and Ir chain atom follow the behavior 
in an infinite monowire very closely as seen in Fig. [Sja) 
and (b). 

The Pt trimer shows a trend similar to Ir upon stretch- 
ing, only with smaller moments — this is in contrast to an 
infinite Pt MW, where the spin moment is zero for a large 
interval of interatomic distances without SOC. It is resur- 
rected upon including SOC, but only when the magneti- 
zation points along the x-axis, an effect coined as colossal 
magnetic anisotropyj^ In contrast, for the suspended Pt 
trimer, the moments of the trimer atoms are non-zero al- 
ready without SOC for a wide range of dc- Upon includ- 
ing SOC the spin moments do not change dramatically 
and the moments are non-zero for both magnetization 
directions in the whole range of AL, which manifests the 
subtle magnetism of Pt chains and its sensitivity to local 
environment. 

In the critical case of Pt, we analyze in more detail the 
influence of the contact geometry on the spin moments 
of the trimer. At the stretching of AL = 1.8 a.u. we per- 
formed a calculation in which we artificially quenched the 
spin moments in the contacts by applying a small mag- 
netic field inside the muffin-tin spheres of the Pt atoms, 
situation referred to as "B2" in the following (see Fig. 2). 
We observe that /i^ drops significantly from 0.55/iB to 
0.18/i_B, while fig drops from 0.89/i_B to 0.54^^ (Fig. 3(c)) 
— a value, very close to that reported in Ref. [1J|, where 
contacts were modelled by infinite bulk electrodes. How- 
ever, quenching the spin moment in the deeper parts of 
the contacts only ("Bl"-case, Fig. 2) — the situation 
which is probably closer to that in real experiments — 
almost does not affect the values of fig and fig. 

Finally, we compare the results obtained within our ID 
geometry of the electrodes with the situation of surface- 
like electrodes, as shown in Fig. 2(a). For a separation 
of AL — 2.4 a.u. fig is reduced in the 2D geometry by 
only 25% from l.OSfiB in the wire geometry to 0.79/^s, 
while the moment of the apex atom is affected stronger, 
dropping from 0.64/l(_b to 0.26/^b (Fig- 3(c)). In a more 
realistic situation the coordination of the contact atoms 
next to the apex atoms is more reduced and the moments 
of the trimer 's atoms will increase approaching the values 
obtained in our wire geometry. However, we have to con- 
clude that magnetism in suspended short Pt chains with 
their small moments is extremely sensitive to slightest 
changes in the contact geometry. Therefore, observing 
this magnetism experimentally will be a hard task, even 
considering the cluster's large magnetic anisotropy ener- 
gies. 



(a) 



> 

B 

< 



(b) 

< 





1 


1 
















-♦ 


-4 


K» = 




— • 






- 

y 


1 


1 


O— O c-atom 
•— • a-atom 
--■ MW 
•— • trimer 







-3' 
-60 
-90 
120 



1 1 1^ 1 1 1- 










1 


1 D-^ 1 


1 



(c) 







> 

g 

< 



1 1 1 

^ 

• — 1 




2D 






- ^x^^ 

~_ Pt ^ 

1 


1 ^ 


S <> 

\ — <> 



AL (a.u.) 



FIG. 4: (color online) Magnetic anisotropy energies of sus- 
pended trimers of (a) W, (b) Ir and (c) Pt as a function of the 
stretching of the contacts AL. The total trimer's anisotropy 
Ai?*"^' (black half-filled circles) is decomposed into the con- 
tribution from the central A_E'' (red open circles) and apex 
atoms Ai?" (blue filled circles) (for definitions of the quanti- 
ties see text). Black dashed line stands for the MAE of the 
infinite MWs (per atom). For Pt squares stand for the AiJ" 
(filled) and AiS"^ (open) tetramer values, while triangles up 
and triangles down stand for the corresponding values in B2 
and 2D cases. For details see text. 



IV. MAGNETIC ANISOTROPY ENERGIES 

Recent theoretical studieSfiii^ which are in accor- 
dance with experiments clearly state that the in- 
teratomic distance dc in suspended chains of hd ele- 
ments consisting of several atoms can reach large val- 
ues of 5.0—5.5 a.u. upon stretching, which corresponds 
to AL w 2 — 3 a.u., c.f. Fig.[TJ Therefore, in a real Pt, Ir, 
or W break junction experiment the central atom will re- 
veal strong fingerprints of spin-polarization at sufficiently 
small temperatures. In order to study the thermal sta- 
bility of magnetism in the suspended trimers, next we 
concentrate on their magnetic anisotropy energies. 

In a transport break junction experiment, the data is 
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obtained by averaging over thousands of measurements, 
which differ from each other by the thinning history of 
the wire and the contact geometry. Therefore, in order 
to analyze this data with respect to the magnetism of the 
suspended chain, it is necessary to disentangle the contri- 
butions of intrinsic trimer's magnetic anisotropy energy 
Ai5*''' from that originating from the contacts. In our cal- 
culations, we do this by switching SOC off in the contact's 
atomic spheres. Moreover, we individually switch SOC 
on and off in the apex and central atoms, to determine 
their separate contributions, AE"^ and AE'^, respectively, 
to the total MAE of the trimer, A£;'". For W with essen- 
tially one magnetic atom in the trimer we define AE^^^ as 
AE''+2AE'', while for Ir and Pt we define the total MAE 
of the trimer per atom as AE*'" = {AE" + 2AE°')/3^ 

The absolute value of trimer's MAE, shown in Fig. 4 
with a black solid line for W, Ir and Pt, reach large ab- 
solute values of 10 to 30 meV per atom at most electrode 
separations. These values, which are one to two orders of 
magnitude larger than those in most of transition-metal 
nanostructures, give us confidence that the magnetism of 
suspended trimers can be tackled experimentally. For all 
three elements the calculated values of anisotropy ener- 
gies AE*"" are of the order of magnitude of those pre- 
dicted for infinite monowires (dashed line in Fig. 4 and 
Refs. [l|Tl|). Interestingly, the trimer's MAE displays a 
non-trivial dependence on AL quite different from that 
in ideal chains for Ir and W, while for Pt both anisotropy 
energies are close in their trend. 

For W, Fig. 4(a), the only contribution to Ai?*''' stems 
from the strongly spin-polarized central atom. Start- 
ing at about 35 meV and a magnetization along the 
chain, there is a switch to a perpendicular direction upon 
stretching by AL = 0.6 a.u. Upon further increasing the 
contact separation, the perpendicular magnetization is 
stabilized by a sizable MAE of 10—20 meV. A switch of 
the magnetization direction can be also observed for Ir at 
significant stretching of around AL = 3.0 a.u., Fig. 4(b). 
In general, for Ir, the behavior of Ai?*''' upon stretching 
is quite smooth, and the ^-direction of the magnetiza- 
tion is favored by considerable 15—30 meV over a wide 
interval of contact separations. 

Detecting traces of magnetism in Pt break junctions at 
small stretching will present a significant challenge as can 
be seen from Fig. 4(c). In the regime of AL < 1.3 a.u. the 
value of AE^" is on the order of a few meVs, posing the 
question of whether the sensitive magnetization will sur- 
vive temperature fluctuations. Moreover, in this regime, 
as pointed out above, the spin moments of the apex and 
central atoms depend strongly on contact geometry and 
details of the thinning process. This will influence the 
contribution of these atoms to Ai?*''', causing frequent 
measurement-to-measurcmcnt oscillations in its sign and 
magnitude. Beyond a distance of AL — 1.5 a.u. the well- 
developed magnetization of the trimer is pointing rigidly 
along the chain axis with a MAE of about 15 meV — a 
value somewhat smaller than that of an infinite MW of 
corresponding interatomic distance, c.f. Fig. 4(c). 



As a general trend, we observe in Fig. 4 competing con- 
tributions to the MAE from central and apex atoms for Ir 
and Pt atomic junctions. At most separations a positive 
value from the apex atom, AE"" > 0, favors a perpendic- 
ular magnetization direction, while negative contribution 
from the central atom, AE"^ < 0, forces the magnetiza- 
tion along the trimer's axis. Due to this competition the 
trimer's MAE behaves qualitatively differently from that 
in the infinite MW as a function of the interatomic dis- 
tance for Ir. It significantly quenches the central atom's 
MAE, so that the resulting total anisotropy is on average 
smaller than anticipated from the infinite MW both for 
Ir and Pt. 

The origin of this competition is the different local 
symmetry of apex and central atoms and it exists also 
in longer chains or chains made of other elements. To 
demonstrate this point we additionally calculated the 
MAE of a suspended four-atom chain of Pt atoms with 
the contacts as in Fig. 1. The AE°' and AE'^ anisotropy 
energies for the two apex and two central atoms, respec- 
tively, are displayed by squares in Fig. 4(c) (per atom) 
and reveal the same effect: AE'^ and AE"^ are close in 
their values but opposite in sign. Due to its large mag- 
nitude, this effect may even compete with exchange in- 
teractions in suspended small chains and lead to a non- 
coUinear magnetic order. With increasing chain length 
the anisotropy energies of the atoms in the center of the 
chain will eventually approach that of an atom in the infi- 
nite monowire (this can be already seen in Fig. 4(c) for Pt 
when going from a trimer to a tetramer), and the whole 
cluster will behave as a superparamagnet. For smaller 
suspended clusters, however, which are much more prob- 
able to occur in an experiment, the effect of competition 
between central and apex atoms will be dramatic. 

We demonstrate that the effect of the competition be- 
tween the apex and central atoms for the MAE value is 
also stable with respect to the geometry of the contacts. 
For this purpose we calculate the MAE of the Pt trimer 
in 2D-configuration, as well as in Bl and B2 cases, and 
present the calculated values in Fig. 4(c). As far as the 
2D-case is concerned, we observe that despite significant 
changes in the spin moments, c.f. Fig. 3(c), the values of 
AE'^ and Ai?^ are very close to those calculated within 
the wire geometry, moreover, they are also opposite in 
sign. Their competition results in rather close values for 
the total AE^" of 7 and 18 meV for the 2D and wire- 
like geometry of the contacts, respectively, both favoring 
the z-direction of the trimer's magnetization. Notably, 
both values are very close to that of 12 meV reported by 
Smogunov et al. in Ref. [l3l for the suspended between 
semi-infinite electrodes Pt trimer, although one has to 
keep in mind the difference between the geometries of all 
three cases. 

In the Bl-situation, in which we do not quench the 
spin moments of the contact atoms of the electrode, the 
MAE is almost not affected compared to the wire setup 
of Fig. 1. On the other hand, in the B2-case, in which we 
quench the spin moments inside the entire electrode, the 
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FIG. 5: (color online) Spin-decomposed local density of states 
(LDOS) for a Pt trimer including SOC for a separation be- 
tween the electrodes of AL = 1.8 a.u. (a) and (b): A3 
{dxz,dyz) and A4 {d^2_y2 ,dxy) 3d-LDOS of the central ("c- 
atom") and apex ("a-atom") atoms for the magnetization 
along the z-axis. (c) LDOS of the central atom for two differ- 
ent magnetization directions {z and r). 



spin moments of the trimer atoms are strongly reduced, 
c.f. Fig. 3(c), and so is the MAE: AE" and AE" drop to 
ahuost half their value. However, the latter two contri- 
butions to the total MAE AE*^" are still of opposite sign, 
proving the generality of the phenomenon of competing 
anisotropics for this type of atomic junctions. The latter 
effect is responsible for the value of the total MAE AE^" 
of 9 meV, smaller than the value of 20 meV obtained 
within the wire-geometry. 

To understand the effect, we analyze the local density 
of states (LDOS) and choose the Pt trimer as an example. 
In an infinite MW, the position of the localized d-states 
with respect to the Fermi energy (Ep) is responsible for 
the formation of the orbital moment and the direction 
of the magnetization ] -'^ 1 ,28,29 According to the symmetry 



these d-states can be subdivided into A3 {dxz,dyz) and 
A4 {dr^2_y2,dxy) groups, not taking into account Ai sd- 
hybridiii In a Pt (Ir) infinite MW, the position of the X- 
edge (F-edge) of the A4 degenerate dxy and dx2-y2 bands 
at Ep is responsible for the easy magnetization direc- 
tion along the chain i^^'^^ In a trimer at larger stretching, 
the 5d-LD0S of a central atom qualitatively follows that 
of an atom in an infinite monowire (not shown), which 
explains the same sign of the MAE (Fig. |4]). A small 
splitting between the d^^^-y^ and dxy states due to the 
presence of the leads can be already seen for the cen- 
tral atom (Fig. 5(a)), however, for the apex atom the 
effect is dramatic and has crucial consequences for the 
apex MAE. Interaction with the contacts locally breaks 
the Cioo symmetry of an infinite chain which results in 
large splitting between the dxy- and d2;2_j,2-orbitals in 
the A4-band and their shift towards the lower energies 
(Fig. 5(b)). On the other hand, the degeneracy between 
the dxz and d^z-states (As-band) is locked by symmetry. 
This leads to the rearrangement of the states around Ep 
and as a result the in-plane direction of the apex spin 
moment becomes favorable (Fig. 4). 

In general, it should be possible to deduce the pre- 
dicted switches in the magnetization direction of the sus- 
pended chains from transport measurements in these sys- 
tems. As an example, in Fig. 5(c) we present the LDOS 
of the central atom in a Pt trimer at AL — 1.8 a.u. Upon 
changing the magnetization direction significant changes 
in the LDOS around the Fermi energy Ep can be ob- 
served (shaded area in Fig. 5(c)), which will inevitably 
result in large changes of the experimentally measured 
conductance. As far as the 5d TMs are concerned, gi- 
ant values of the ballistic anisotropic magnetoresistance 
should be achievable due to strong spin-orbit coupling 
in these metals, which is able to modify the electronic 
structure significantly in response to the changes of the 
magnetization direction in real space. 



V. CONCLUSIONS 

By performing ab initio calculations of suspended 
trimers of W, Ir and Pt including structural relaxations 
as a function of the separation between the leads, we 
demonstrate that the chain atoms develop significant spin 
moments upon stretching. We investigate the stabil- 
ity of these spin moments via evaluating the magnetic 
anisotropy energy of the trimers. Our calculations re- 
veal large MAE of the whole trimers of the order of 10 
to 30 meV per atom. Interestingly, we observe that the 
total MAE presents a competition between large contri- 
butions from the apex and central atoms. We argue that 
this effect is general and can occur in suspended chains 
of different elements and different length, leading to non- 
trivial real space textures of magnetic anisotropy energy 
which might even lead to complex magnetic ordering in 
atomic-scale junctions. 

Under the condition of large predicted values of MAE 
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huge magnetic fields would be required to change the 
magnetization of the chains, therefore, new ways to 
achieve this goal have to be established. One of the pos- 
sibilities lies in using the ability to control the magneti- 
zation direction in the junction by changing the distance 
between the leads. We show that such switches of the 
magnetization can happen and will result in strong fea- 



tures in the measured conductance. Distinguishing these 
features from the changes in conductance due to chain 
elongation or atomic rearrangements presents a consid- 
erable challenge. 
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